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Abstract 



We investigate the dynamical processes taking place in nanodevices driven 
by high-frequency electromagnetic fields. We want to elucidate the role of dif- 
ferent mechanisms that could lead to loss of quantum coherence. Our results 
show how the dephasing effects of disorder that destroy after some periods 
coherent oscillations, such as Rabi oscillations, can be overestimated if we 
do not consider the electron-electron interactions that can reduce dramati- 
cally the decoherence effects of the structural imperfections. Experimental 
conditions for the observation of the predicted effects are discussed. 
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Recent advances in laser technology make possible to drive semiconductor nanostruc- 
tures with intense coherent ac-dc fields. This opens new research fields in time-dependent 
transport in mesoscopic systems and puts forward the basis for a new generation of ultra- 
high speed devices. Following these results several works have been devoted to the analysis 
of the effects of time-dependent fields on the transport properties of resonant heterostruc- 
tures |2|-§J, and to exploit the application in ultrafast optical technology: high-speed optical 
switches, coherent control of excitons, etc. f>]||. 

Within this context, we are interested in the decoherence processes producing the ob- 
served fast dephasing of coherence phenomena in semiconductor superlattices (SL's) and 
more specifically in the interplay between the growth imperfections (disorder) and many- 
body effects as the electron-electron (e-e) interaction. The interplay between the effects 
of disorder and many-body effects on electronic properties is a long-standing problem in 
solid-state physics. Probably one of the most promising way to gain insight into this in- 
tricate problem is to combine the actual state-of-the-art of the Molecular Beam Epitaxy 
(MBE), which allow us to grew samples with monolayer perfection and consequently with 
well-characterized disorder, with coherent oscillations that are extremely sensitive to imper- 
fections and nonlinear effects. 

The oscillations of a two level system between the ground and excited states in the 
presence of a strong resonant driving field, often called transient nutation or Rabi oscillation 
(RO), are discussed in textbooks as a topic of time-dependent perturbation theory. Consider 
a two state system with ground state energy E and excited state Ei in the presence of a 
harmonic perturbation. If the frequency of the perturbation matches roughly the spacing 
between the two levels, the system undergoes oscillations with a frequency Qr which is much 
smaller than the excitation frequency uj ac . This Rabi frequency depends on the mismatch 
5u) = {Ei — E )/h — uj ac between the level spacing and the excitation frequency, and on 
the matrix element Eiq of the perturbation Qr = (j)u 2 + |i<io| 2 /^ 2 ) ^ . If we start with the 
system initially in the ground state, transitions between the ground and the excited state 
will occur with a period Tr = 2'k/VIr as time evolves. 
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Semiconductor SL's present Bloch minibands with several states each one; thus it is not 
clear whether they can be correctly described as a a pure two state system. We should 
also take into account the presence of imperfections introduced during growth processes 
and scattering mechanism as e-e interactions. Interface roughness appearing during growth 
in actual SL's depends critically on the growth conditions. For instance, protrusions of 
one semiconductor into the other cause in-plane disorder and break translational invariance 
parallel to the layers. To describe local excess or defect of monolayers, we allow the quantum 
well widths to fluctuate uniformly around the nominal values; this can be seen as substituting 
the nominal sharp width by an average along the parallel plane of the interface imperfections. 
Our approximation is valid whenever the mean-free-path of electrons is much smaller than 
the in-plane average size of protrusions as electrons only see micro-quantum- wells with small 
area and uniform thickness. Therefore, in the following we will take the width of the nth 
quantum well to be a(l + We n ), and the width of the nth barrier as 6(1 + We n ) where W is a 
positive parameter measuring the maximum fluctuation, e n 's are distributed according to a 
uniform probability distribution, V(e n ) = 1 if |e n | < 1/2 and zero otherwise, a is the nominal 
quantum well width and b is the nominal quantum barrier width. Even with its rather 
satisfactory degree of success, many-body calculations have difficulties that, in some cases, 
may complicate the interpretation of the underlying physical processes. Presilla et al. @ 
suggested a new treatment of these interactions that, loosely speaking, could be regarded as 
similar to Hartree-Fock and other self-consistent techniques, which substitute many-body 
interactions by a nonlinear effective potential. More recently || we have proposed a new 
model where the nonlinear interaction is driven by a local field instead of the mean-field 
approach used by Presilla and co-workers. 0. 

Here we present a new model that solves self-consistently the Poisson and the time- 
dependent Schodinger equations, and we compare the results with a local treatment based 
in the non-linear Schrodinger equation. 

The band structure and the potential at flat band is computed by using a finite-element 
method. The eigenstate j of the band i with eigenenergy is denoted as ipf\x). A good 

3 



choice for the initial wave packet is provided by using a linear combination of the eigenstates 
belonging to the first miniband. For the sake of clarity we have selected as the initial wave 
packet ty(x,0) = ip^(x), although we have checked that this assumption can be dropped 
without changing our conclusions. The subsequent time evolution of the wave packet ^f?(x, t) 
is calculated numerically by means of an implicit integration schema designed for consider 
time-dependent fields [f|. 

The envelope-functions for the electron wavepacket satisfies the following quantum evo- 
lution equation 



.+ d$(x,t) 



n 2 d 2 Tr , ; 

+ V NL (x,t) 



2m* dx 2 



*(x,t), (1) 



where x is the coordinate in the growth direction of the SL. We consider two approaches 
to the nonlinear potential Vnl{x, t) in Eq.( [1]). On the one hand we take into account the 
model described, in other context, in || for our problem, where Vnl(x, t) is 

V NL (x,t) = V(x) -eF Ac xsm(u Ac t) + a loc \^(x,t)\ 2 , (2) 

and V(x) is the potential at flat-band, F A c and ujac are the strength and the frequency of 
the ac field respectively, and all the nonlinear physics is contained in the coefficient ai oc , 
which we discuss below. There are several factors that configure the medium nonlinear 
response to the tunneling electron. We want to consider only the repulsive electron-electron 
Coulomb interactions, which should enter the effective potential with a positive nonlinearity, 
i.e., the energy is increased by local charge accumulations, leading to a positive sign for a loc . 

On the other hand, we have considered a different approach by solving self-consistently 
the Schrodinger and Poisson equations obtaining a Hartree-like potential. In this context, 
the non-linear potential is, 

V NL (x, t) = V(x) - eF Ac x wa.(u AC t) + a se ifV H (x, t), (3) 

where now Vh it is obtained by solving the Poisson equation for the density of charge 
\*&(x, t)\ 2 , and a se if is the coupling parameter. 



We present here results for a SL with 10 periods of 100 A GaAs and 50 A Gao.7Alo.3As 
with conduction-band offset 300 meV and m* = 0.067m, m being the free electron mass. 
To illustrate the effects of the nonlinear coupling we show in Fig. [I] the conduction-band 
profile for a perfect SL (W — 0) at t — 0.4 (lower) and 1.2 ps (upper) when the ac field is 
tuned to the resonant frequency uj ac = co> res ~ 24 THz, for (a) the linear case and modeling 
the e-e interaction with (b) the self-consistent method (a se ;/ = 10~ 3 ) and with (c) the 
local model («/ oc = 10). To show the effects of the interface roughness we plot in Fig. |2| 
the probability of finding an electron, initially situated in ip^ (x) , in the state ip[ 5 \x) as a 
function of time when the ac field is tuned to the resonant frequency u; ac = u) rcs ~ 24 THz 
with (a) W = (perfect SL) and (b) W = 0.03 (imperfections around one monolayer). 
These results suggest the existence of a characteristic scattering time Tdi S related to the 
amount of disorder in the sample, after which RO's are destroyed by disorder. In Fig. ||] we 
plot the probability of finding an electron, initially situated in ipQ 5 \x), in the state ip[ 5 \x) as 
a function of time when the ac field is tuned to the resonant frequency, for different values 
of the nonlinearity coupling (a) a se if = 5 x 10~ 5 , (b) 10~ 4 , (c) 5 x 10~ 4 and (d) 10~ 3 . The 
results for the local model are very similar. When we compare this picture with Fig. ^|we see 
the process of vanishing of the RO's are completely different. In the second case, the effects 
are the same for any time, then we could not speak about a dephasing time, apparently 
we only modified the electronic structure and then we are decreasing the resonant coupling 
between the external ac-field and the Bloch bands. Nevertheless the main goal of this work 
is to show how the nonlinearity effects can reduce the dephasing effects introduced by the 
growth imperfections. In Fig. [| we plot the occupation probability of the state ifj[ 5 \x) as 
a function of time considering imperfections about one monolayer (W = 0.03) for (a) the 
linear case and considering together with the imperfections the e-e interaction with the self- 
consistent model (b) a se if = 10~ 4 and with the local one (c) a\ oc = 5. We can see clearly 
how nonlinearity prevents the dephasing effects introduced by the imperfections allowing 
the observation of Rabi oscillations during larger coherence times. These theoretical results 
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are completely consistent with recent experiments in transport properties of intentional 



disordered superlattices with doped and undoped superlattices p3| , |TT]| , where is showed 
than the Coulomb interactions could be the responsible of the suppression of disorder effects 
leading to quasimetallic behavior at low temperatures when the doping of the samples is 
increased. 

In summary, we have shown how the dephasing effects of disorder are dramatically re- 
duced when we consider the e-e interaction. We have studied two different models to in- 
troduce the non-linear interaction and the results are very similar. Our results shows that 
it is possible to enlarge the dephasing times and, consequently the number of periods of 
coherence oscillations of electrons in SL's. In semiconductor heterostructures this can be 
done by increasing the doping or with very intense laser excitation fields. It goes without 
saying that to develop new devices for THz science it is crucial to understand how to control 
and enlarged the coherence times. We think that the nonlinear effects could be the key to 
solve this problem. Further work along these lines is currently in progress. 
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FIGURES 



FIG. 1. Conduction-band profile. 

FIG. 2. The probability of finding the electron in the state ip^\x) as a function of time for 
different amounts of disorder. 

(5) 

FIG. 3. The probability of finding the electron in the state (x) as a function of time for 
different values of the nonlinear coupling. 

FIG. 4. The probability of finding the electron in the state tp[ 5 \x) as a function of time 
considering the interplay between imperfections and nonlinearities. 
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